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Abstract 

We perform an update of the next-to-leading order calculation of the rate for Higgs boson 
production in association with two jets. Our new calculation incorporates the full analytic result 
for the one-loop virtual amplitude. This new theoretical information allows us to construct a code 
including the decay of the Higgs boson without incurring a prohibitive penalty in computer running 
time. Results are presented for the Tevatron, where implications for the Higgs search are sketched, 
and also for a range of scenarios at the LHC. 
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I. INTRODUCTION 



In the coming years the hadron colliders at Fermilab and CERN will focus on the hunt 
for the Higgs boson. The large data sample currently being collected at the Tevatron will 
certainly lead to improved limits on the Higgs mass [1], or even evidence for its existence. The 
Large Hadron Collider has the potential to confirm the existence of the Higgs boson {2, ^] 
between the lower .hnit set by LEP Q and the npper bonnd snggested by pertnrbative 
unitarity |5|, |6j. 

Such claims are based on detailed analyses that clearly require reliable theoretical pre- 
dictions for the production cross sections and characteristics. It is well-known that leading 
order predictions for such quantities, based on tree-level Feynman diagrams alone, are not 
sufficiently trustworthy for this purpose. The calculations are plagued by large uncertainties 
in their overall normalization and moreover, important kinematic effects are often missed. 

In this paper we present results for the production of a Higgs boson in association with 
two jets. Our calculation is performed at next-to-leading order (NLO) using an effective 
Lagrangian to express the coupling of gluons to the Higgs field Q], 

C^ = ^HtiG^G^. (1) 

where the trace is over the color degrees of freedom. At the order required in this paper, 



the coefficient C is given in the MS scheme by 



C=£(l + % 8 )+0{a? 8 ). (2) 
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Here v is the vacuum expectation value of the Higgs field, v = 246 GeV. 

This Lagrangian replaces the full one-loop coupling of the Higgs boson to the gluons via 
an intermediate top quark loop by an effective local operator. The effective Lagrangian 
approximation is valid in the limit m# < 2m t and, in the presence of additional jets, 
when the transverse momenta of the jets is not much larger than the top mass m t 10] . A 
commonly used improvement of the effective Lagrangian approximation is to multiply the 
resulting differential jet cross section by a ratio R given by, 

D °"finite m t (99 -> H) 

K = 1 ruv > ' 6 ) 

cr mt ^oo{gg -»■ H) 

where a(gg —> H) is the total cross section. Setting x = Am1/m 2 H the correction for the 
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finite mass of the top quark in the region x > 1 is 8], 



R 



3x ( , . 



-I 2 



sin 



(4) 



This rescaling is known to be an excellent approximation for the Higgs + 2 jet rate, see 



Rcf. 



10] and references therein. However for the case of Higgs + 1 jet it has been found 



that the effect of bottom quark loops and additional electroweak diagrams can also be 
important 11] and these effects should also be included. Our numerical results for the Higgs 
cross section will not include the rescaling of Eqs. ( I3l4l) . 



II. NEW FEATURES OF THIS PAPER 



The phenomenology of the production of a Higgs boson in association with two jets has 
been presented in Ref. [l2j for the case of the LHC operating at y/s = 14 TeV. The NLO 
analysis in that paper was based on real matrix elements for the Higgs+5 parton processes 



given in Ref. 



13] , supplemented by the results of Ref. [l4j, |l5| in the cases where these latter 



results lead to more efficient code. In Ref. [12] the virtual matrix element corrections for the 
Higgs + 4 parton process were taken from Ref. 



161 ] . For the Hgggg and Hqqgg sub-processes 
the virtual corrections were based on a semi-numerical technique 17| . whilst the matrix 
elements squared for the one-loop processes Hqqq'q' and Hqqqq were given analytically in 
Ref. [16]. 

In the three years since Ref. [3] was published a great deal of effort has been devoted 
to the analytic calculation of one-loop corrections to Higgs + n-parton amplitudes, with 
particular emphasis on the n = 4 amplitudes which are relevant for this study. The complete 
set of one-loop amplitudes for all Higgs + 4 parton processes is now available and analytic 
expressions can be found in the following references: 



Hgggg: Refs. [lS^; 



Hqqgg: Refs. [23j,|24j; 



HqqQQ: Ref. [23] 



These new analytic results have now been included in the MCFM package, version 5.7 
(which may be downloaded from mcfm.fnal.gov), leading to a considerable improvement 



in the speed of the code. For the processes involving two quark-antiquark pairs, the matrix 
elements sqnared given in Ref. fl are implemented in MCFM, rather than the amplitndes 
of Ref. [23J , because they lead to faster code. The values of the amplitudes calculated by the 



new analytic code and the previous semi-numerical code [12|] are in full numerical agreement 
for all amplitudes. 

The improvement in the performance of our numerical code means that it is appropriate 
to revisit the phenomenology of Higgs + 2 jet production and to extend it in a number of 
ways. The improvement in the speed of the code means that it is possible to include the 
decays of the Higgs boson, specifically for the processes: 



h x + h 2 -» H + j± + j 2 r+ + r + j x + j 2 

h 1 + h 2 ^ H + j 1 +j 2 ^b + b + j l + j 2 

h 1 + h 2 ^H + Jl + j 2 W~ + W + + h + j 2 

U- v + e + 
-»• e~ + v 



(5) 
(6) 

(7) 



h x + h 2 H + ji + j 2 Z + Z + 3l + j 2 

U e~ + e+ (8) 
-> IT + /i + 

where /ii, /12 represent partons inside the incident hadron beams. All four of these processes 
are included in MCFM v5.7. 



III. PARAMETERS 



Throughout this paper we make use of the MSTW2008 parton distribution functions [25] , 
using the LO fit (a s (Mz) = 0.13939 and 1-loop running) for the lowest order calculation 
and the NLO fit (a s (M z ) = 0.12018 and 2-loop running) at NLO. The W mass and width 
are chosen to be, 

m w = 80.398 GeV, T w = 2.1054 GeV . (9) 



The mass is taken from Ref. The total width given in Eq. Q is derived from the 



measured branching ratio for W —> £u, 10.80 ±0.09% [26( by using a lowest order calculation 
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°lo [fb] 


&NLO [fb] 


°lo [fb] 


<^ NLO [fb] 


Hi; 


Igs + jets 


1.25 


1.98 


1.25 


2.05 


Hi; 


Igs + 1 jets 


0.84 


1.16 


0.74 


1.07 


Hig; 


t,s + > 2 jets 


0.35 


0.48 


0.28 


0.39 



TABLE I: Cross section for Higgs + jet production and decay into W~{— > fi~P)W + (— > ve + ) at 
yfs = 1.96 TeV for Mh = [i= 160 GeV. In the second and third columns, only the cuts of Eq. (jlip 
are applied. For the results in the final two columns the more stringent cut, |?7jet| < 2 is applied, 



in order to allow a comparison with Ref. 
of the partial width, 



29]. 



r(w -> to) 



(10) 



y/2 6vr 

This ensures that our calculation incorporates the best possible value for the W branching 
ratio which is determined to about 1%. The values of the total Higgs width are taken from 
the program hdecay 27J, version 3.51. 



To define the jets we perform clustering according to the algorithm 
definitions detailed further below. 



28 



with jet 



IV. TEVATRON RESULTS 



We use a very simple set of inclusive cuts, with no requirements on the Higgs boson decay 
products, 

(jet) > 15 GeV, |^ jet | < 2.5, i? je t, je t > 0.4 . (11) 

At the Tevatron the search for the Higgs boson has been divided into jet bins. To set 
the stage for this we show in Table [I] the expected cross section in each bin due to the gluon 
fusion mechanism. The parameter jj, is the renormalization and factorization scale, which 
we set equal to tuh here. We note that next-to-next-to-leading order (NNLO) results for the 
Higgs + jet cross section are given in [291], based on the earlier calculations in Refs. 30-32]. 
From table [J columns 3 and 5, we see that the Higgs + > 2 jets bin constitutes about 13% 
of the cross section for \r)j et \ < 2.5 and 11% with \rjj e t\ < 2. 

It is interesting to compare the number for the fraction of Higgs + > 2 jet events (|%t| < 
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m H [GeV] 


150 


160 


165 


170 


180 


T H [GeV] 


0.0174 


0.0826 


0.243 


0.376 


0.629 


<?lo [fb] 

&NLO [fb] 


u - 3zy -45% 


0.476+^% 


n qqi +92% 
0.458+36% 


Qnci+92% 
0.422+S 


245 +91% 
n ^4 c i+ 37 % 

U.J4D_ 31% 


Finite correction, R 


1.098 ± 0.003 


1.113 ±0.003 


1.122 ±0.004 


1.130 ±0.005 


1.149 ±0.005 



TABLE II: Cross section for Higgs + 2 jet production and decay into W~{— > fi~u)W + (— > ve + ) at 
yfs = 1.96 TeV. Only the cuts of Eq. ([lip are applied. The correction factor for each Higgs mass, 
given by Eq. is also shown. 



2) with the percentage extracted from Table 2 of [291 ]. which is quoted as 4.9%. Our number 
is deficient in that it does not include NNLO corrections to the Higgs + jet rate. Our 
calculation treats all jet bins consistently at NLO. The inclusion of the NNLO correction to 



the Higgs + jet bin will reduce our number. On the other hand, the calculation of Ref. [29 1 
is deficient because it does not treat all bins consistently at NNLO, i.e. it does not include 
NNLO corrections for the Higgs ±1 jet rate or NLO±NNLO effects for the Higgs + > 2 
jet rate. We roughly estimate that including the NNLO effects in the Higgs ±0 jet bin 
would move our central value from 11% to 10%. Overall, because the corrections are quite 
substantial, the theoretical estimate of the fraction of events in the Higgs + > 2 jet bin is 
quite uncertain. 

Despite the fact that the fraction of events in the Higgs + > 2 jet bin is small, it is 
important because the associated uncertainty is large. We investigate this issue in Table HT1 
where we give the cross section for the process of Eq. (]7|) using a selection of values for the 
Higgs mass of current interest for the Tevatron. In the table we give the results for the leading 
order and next-to-leading order cross sections, calculated using LO and NLO MSTW2008 
PDFs respectively. For the range of Higgs masses considered, the QCD corrections increase 
the cross section by approximately 40% (for the central value, \x = m#). The theoretical 
error is estimated by varying the common renormalization and factorization scale in the 
range, mn/2 < fi < 2m#. As can be seen from the table, even though including the next- 
to-leading order corrections leads to a considerable improvement in the theoretical error, the 
remaining error is still quite sizeable. We do not include a factor to correct for the finite top 
mass, but in order to facilitate comparison with other calculations we also tabulate this factor 
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R (computed using Eq. (ffl ) using a value for the top quark mass of m t = 172.5 ± 2.5 GeV. 

In the spirit of Ref. [29], we can now estimate the theoretical uncertainty on the number 
of Higgs signal events originating from gluon fusion. By using the fractions of the Higgs 
cross section in the different multiplicity bins taken from Ref. [33], we can update Eq. (4.3) 



of Ref. 



29| (for a Higgs boson of mass 160 GeV) with, 



A7V T al(scale) = 60% ■ CS) + 29% ■ KS) + 11% • (US) = KKS (12) 

-< > signal 

Only the uncertainty on the Higgs + > 2 jet bin has been modified, using the results from 
Table [TTJ The corresponding determination using the LO uncertainty in the Higgs + > 2 



jet bin is (+20, 0%, —16.9%) 29], so this represents a modest improvement in the overall 
theoretical error. 

The correspondence of our results with those of Anastasiou et al. is somewhat obscured 
by the fact that the total Higgs width used in Ref. [29| is about 7% smaller at tuh = 160 GeV 
than the value given in our Table HIl Taking this fact into account and including the finite 
top mass correction tabulated in Table HT1 we find that our NLO Higgs + 1 jet and LO Higgs 
+ 2 jet cross sections in Table [J are in agreement with the corresponding numbers (1.280 
and 0.336 fb) from Table 2 of Ref. [2J 



A. Effect of additional search cuts 

We also investigate the behaviour of the LO and NLO predictions in the kinematic region 
relevant for the latest Tevatron Higgs exclusion limits. Therefore, in addition to the jet cuts 
above, we also consider cuts on the decay products of the W/W* that are produced by the 



Higgs boson. These cuts correspond very closely to a recent CDF analysis 34j, although 
the treatment of lepton acceptance is simplified. 

• One of the leptons from the W decays (the "trigger" lepton, £i) is required to be 
relatively hard and central, p^ 1 > 20 GeV, I?/ 1 ! < 0.8 whilst the other (£2) may be 
either softer or produced at slightly higher pseudorapidity, p e t 2 > 10 GeV, |?/ 2 | < 1.1. 

• The invariant mass of the lepton pair is bounded from below (to eliminate virtual 
photon contributions), m^ 2 > 16 GeV. 
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Each lepton must be isolated. Any jet found by the algorithm that lies within a rj — <p 
distance of 0.4 from a lepton should have a transverse momentum less than 10% of 
that of the lepton itself. 

The missing transverse momentum - in our parton level study, the sum of the two 



neutrino momenta - is constrained using the lf t spec variable defined by 34 1 



$ spcc = # sin [min (A0, |)] . (13) 



A(j) is the distance between the 1ft vector and the nearest lepton or jet. We require 
that $ t spec > 25 GeV. 

In Figure [TJ we see the scale dependence of the LO and NLO cross sections for ran = 
160 GeV. The upper two curves show the case of the minimal set of cuts in Eq. (ITT)) and 
the lower curves show the results when including the Higgs search cuts above. Applying 
the additional cuts on the Higgs decay products does not change the scale dependence, 
indicating that the isolation and missing transverse momentum cuts (that are sensitive to 
additional radiation) do not play an important role. Applying the additional search cuts 
does not alter the behaviour of the NLO prediction in the Higgs + > 2 jet bin, so that the 
results presented in the previous section (with no cuts on the Higgs decay products) are 
sufficient to estimate the percentage theoretical uncertainty. 

V. LHC RESULTS 

In order to study the impact of the NLO corrections at the LHC, we adopt a different 
set of cuts to define the jets. The rapidity range of the detectors is expected to be much 
broader, allowing for a larger jet separation too, and we choose a somewhat higher minimum 
transverse momentum, 

p*(jet) > 40 GeV, |r/ jet | < 4.5, i2 jetijet > 0.8 . (14) 

In this section we do not consider the decay of the Higgs boson for the sake of simplicity. 



Since results for this scenario have already been discussed at some length [12] , we restrict 
ourselves to a short survey of the essential elements of the phenomenology at the lower centre- 
of-mass energy, y/s = 10 TeV. We present the scale dependence of the LHC cross section 
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m H [GeV] 


120 


140 


160 


180 


200 


T H [GeV] 


0.0036 


0.0083 


0.0826 


0.629 


1.426 


o"lo [pb] 
onlo [pb] 


1 QQ + 78% 

i-8«_4 % 
±.jo_ 23% 


1 48 +76% 

l.DJ_ 23% 


1 20 +75% 

1 q/^+23% 
l.JD_ 23% 


98 +74% 

1 ic+24% 

i.io_ 23% 


°- 81 -39% 

98 +25% 
u.yo_ 24% 


Finite m< correction, R 


1.060 ± 0.002 


1.084 ± 0.003 


1.113 ±0.004 


1.149 ±0.005 


1.191 ±0.007 



TABLE III: Cross section and uncertainties for Higgs + 2 jet production at y/s = 10 TeV with the 
cuts of Eq. (|14p . The correction factor for each Higgs mass, given by Eq. is also shown. 

for Higgs + 2 jets (m# = 160 GeV) in Figure [2j We have also checked the agreement of our 
calculation with previous results 12| at y/s = 14 TeV, taking into account the different choice 
of parton distribution functions used in that reference. As noted in the earlier paper 12], 
the corrections are quite modest using our central scale choice, /io = /4ff> increasing the cross 
section by approximately 15%. Once again, although the scale dependence is much reduced 
it is still substantial. 

For the sake of illustration we have chosen m# = 160 GeV in the study above. To 
illustrate the effect of the QCD corrections more broadly, in Table II III we give the cross 
sections for Higgs masses in the range 120 GeV < itlh < 200 GeV. It is within this range 
that the Higgs + 2 jet process considered here is of most interest, due to its interplay with the 
electroweak weak boson fusion channel. We observe that the effect of the QCD corrections 
increases from about 5% for m# = 120 GeV to 21% for m# = 200 GeV. Estimating the 
theoretical error in the same way as before, we see that the uncertainty is slightly less at 
the LHC than at the Tevatron. 

It is also interesting to consider the dependence of the cross section on the minimum 
transverse momentum required for the observed jets. Results for several other values of this 
threshold, either side of our default value of 40 GeV, are shown in Table IIVI As can be seen 
from the table, the percentage effect of the NLO corrections on the total rate is practically 
independent of the value of p™ m (jet) in the range studied. 



A. Weak boson fusion 



As noted above, the process studied in this paper produces the same final state as expected 
from Higgs production via weak boson fusion (WBF). Although the electroweak process is 
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FIG. 1: Scale dependence for the Higgs + 2 jet cross section, with the Higgs decay into W~(— > 
(j,~v)W + (— > ve + ), at the Tevatron and using the a central scale (jlq = Mh- Results are shown 
for the minimal set of cuts in Eq. (|lip (upper curves) and for cuts that mimic the latest CDF 
H — > WW* analysis (lower curves). 
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FIG. 2: Scale dependence for the Higgs boson + 2 jet cross section, using the basic set of cuts in 
Eq. and a central scale choice /jlq = run- 
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P r n (jet) [GeV] 


20 
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o"LO [pb] 
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TABLE IV: Cross section for Higgs + 2 jet production at sfs = 10 TeV, with m# = 160 GeV and 
the minimum jet p t allowed to vary from that specified in Eq. (|14|) . 

expected to dominate once appropriate search cuts are employed, the remaining fraction of 
events originating from gluon fusion must be taken into account when considering potential 
measurements of the Higgs coupling to W and Z bosons. 

To address this issue, in this section we present a brief study of the rate of events ex- 
pected using typical weak boson fusion search cuts. In addition to the cuts already imposed 
(Eq. (Ej)), these correspond to, 

l%-%l> 4 - 2 > Vh-Vh<°i ( 15 ) 

where ji and j'2 are the two jets with the highest transverse momenta. These cuts pick 
out the distinctive signature of two hard jets in opposite hemispheres separated by a large 
distance in pseudorapidity. This is illustrated in Fig, [3J, where we compare the distributions 
of the jet pseudorapidity difference (without these cuts) in both gluon fusion and weak 
boson fusion. We note in passing that the shape of this distribution for the weak boson 
fusion process is slightly altered at NLO, whilst the shape of the prediction for the gluon 
fusion process is essentially unchanged. 

In figure H] we show the dependence of the cross section on the c.o.m. energy, from 
= 7 TeV (corresponding to the initial running in 2010-11) to = 14 TeV (design 
expectations). We show the cross section both before and after application of the additional 
weak boson fusion search cuts given in Eq. f)15p . together with the corresponding results 
for the WBF process (also calculated using MCFM [36(). The QCD corrections to both 
processes decrease slightly as \fs is increased, whilst the ratio of the gluon fusion to WBF 
cross sections after the search cuts are applied increases from 20% at 7 TeV to 35% at 
14 TeV. This indicates that, viewed as a background to the weak boson fusion process, the 
hadronic Higgs + 2 jet process is less troublesome at energies below the nominal design 
value. 
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FIG. 3: The jet pseudorapidity difference in gluon fusion (red) and weak boson fusion (blue). The 
NLO predictions are shown as solid histograms, while the dashed lines indicate the LO predictions 
normalized to the corresponding NLO cross sections. 

VI. CONCLUSIONS 

In this paper we have presented phenomenological predictions for the production of a 
Higgs boson and two jets through gluon fusion. These predictions have been made possible 
through the implementation of recent compact analytic results for the relevant 1-loop am- 
plitudes |18l-l24|. The speed with which these amplitudes can be evaluated has enabled us 



12|, with 



to improve upon an existing semi-numerical implementation of the same process 
various decays of the Higgs boson now included. 

We have investigated the behaviour of the NLO cross section at the Tevatron, where 
contributions fr o,„ this channel fornt part of the event sa m p,e tot the .atest Higgs searches Q. 
We find that corrections to the event rate in the Higgs + > 2 jet bin are modest and that 
the estimate of the theoretical error is reduced by approximately a factor of two compared 
to a LO calculation. The resulting error is still rather large, corresponding to approximately 
+40% and —30% across the region of Higgs masses, 150 GeV < m H < 180 GeV. 

For the LHC we have provided a brief study of the behaviour of our predictions for 
collisions at = 10 TeV. We have also performed an analysis of this channel in the 
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context of detecting a Higgs boson via weak boson fusion, where the improved theoretical 
prediction presented in this paper is essential in the long-term for making a measurement 
of the Higgs boson couplings to W and Z bosons. 
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FIG. 4: The yfs dependence of the cross section for mu = 160 GeV at LO (dashed) and NLO 
(solid). Results are shown for the minimal set of cuts in Eq. (I14j) (two upper red curves) and after 
application of the additional WBF Higgs search cuts given in Eq. (I15D (two lower red curves). The 
cross section for the weak boson fusion process is also shown for comparison (four central blue 



curves . 
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